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We examine the Bloch-Peierls-Berry dynamics under a classical nonequilibrium dynamical formulation. 
In this formulation all coordinates in phase space formed by the position and crystal momentum space 
are treated on equal footing. Explicitly demonstrations of the no (naive) Liouville theorem and of the 
validity of Darboux theorem are given. The explicit equilibrium distribution function is obtained. The 
similarities and differences to previous approaches are discussed. Our results confirm the richness of the 
Bloch-Peierls-Berry dynamics. 

1 Introduction 

One of the fundamental dynamical equations in condensed matter physics is the so-called Bloch equation, 
describing one electron moving in a periodic potential in position and crystal momentum space [Q: r = 
^Vk e(k) , Kk = — eV r 4>(r) . Here r S 5R 3 is position, k £ SR 3 crystal momentum, e(k) £ 5? 1 kinetic energy, 
4>(r) £ 5ft 1 the electric potential, e the electric charge of an electron, and h the Planck constant. The index 
I in the gradient operation, V/ = d/dl indicates the corresponding space coordinate. A scalar function 
7~C(r, k), termed the Bloch Hamiltonian, can be defined as 

W(r,k) = l[e(k) + e#r)] . (1) 

In this case, the Bloch dynamics can be rewritten in the canonical Hamiltonian form: r = V^Ti and 
k = —V r TL. The usual incompressible condition, the Liouville theorem, in the six dimensional x r = (r, k) 
phase space follows, 

V x • x = V r • r + V k • k 

= V r • V k H - V k • V r H 

= (2) 

R. Peierls successfully extended the Bloch equations to the case with a weak magnetic field B(r) pQ: 
r = -^Vk e(k), and Kk = — eV r cj)(r) — er x B(r). The Hamiltonian in the form of Eq.lpQ) remains the 
same for this extended dynamics, and the Liouville theorem of Eq.© still holds. Everything appears as 
expected. 

However, the Bloch equations of electron dynamics have been recently modified to include both magnetic 
field B(r) and Berry curvature fi(k) in crystal momentum space 2 . The following extended dynamical 
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equations in the semiclassical limit were obtained, termed Bloch-Peierls-Berry equations in the present 
paper: 

r = iv k e(k) - k x O(k) (3a) 
n 

Hk = -eV r 4>(r) - er x B(r) (3b) 

Interesting, unexpected, and rich behaviors occur in the context of this novel dynamics [21 El- For example, 
it has been found that in six dimensional phase space the current flow becomes compressible, that is, the 
Liouville theorem does not hold Pj: 

V x • x + 0. (4) 

Thus the straight-forward conventional means to Hamiltonian analysis does not appear available. Such a 
situation has been called non-Hamiltonian dynamics P|. 

To remedy this problem of non-Hamiltonian dynamics, Xiao et al. p] introduced a density correction 
factor, denoted J(r, k) = 1 + |B • f2 in the present paper, to force the divergence in phase space to zero: 

V r -(J(r,k)r) + V k -(J(r,k)k) =0 (5) 

Of primary motivation using this technique was in the definition of a phase distribution function consistent 
with the Liouville theorem, because the Liouville theorem states that phase volume in a canonical Hamil- 
tonian formulation is conserved along it's trajectory. This has some good properties. For example, with 
the Liouville theorem a distribution function over a system phase, x, will be time-invariant. To verify this 
property we begin with the total time derivative of a distribution /(x, t) in phase space, ^ = + x • V x / 
and the continuity equation, + V x ■ (x/) = 0. When the Liouville theorem holds, V x • x = 0, the 
continuity equation is identical to the time-invariant distribution, df jdt = 0. 

Disagreement on the treatment in Ref.p] exist in literature 01 03 H3 Ej • For example, in a comment 
from Duval et al. , p] , an alternate approach was pointed out such that a local canonical form, and likewise 
a Hamiltonian description, is achievable. In addition to the question of mathematical formulation, there 
is a real issue of possibly different physical consequences in different approaches [S]. 

In the present paper we study the problem from a completely different point of view to see how the 
Hamiltonian like structure and the equilibrium distribution function emerge from a nonequilibirum process: 
the Darwinian dynamics [Dj- We will show that a Hamiltonian, or energy function, naturally emerges. In 
our demonstration it is clear that the Berry phases due to the magnetic field and the Berry curvature of 
the momentum space appears in equal footing. Furthermore, such a procedure provides a straightforward 
discussion on the equilibrium distribution: a nonequilibrium setting provides a natural way to select the 
steady state distribution. Not all of our results are new. Nevertheless, our demonstration appears to provide 
a clear, consistent, but completely classical starting point to detail the system characteristics as previously 
described while possessing no proclivity toward an incompressible or canonical phase. For example, a 
direct solution to the Fokker-Planck equation has been detailed. This solution is found independent of the 
Liouville theorem and may offer a unique insight into a probability distribution. 

We begin the remainder of the paper in section 2 by first describing the generic decomposition of 
Darwinian dynamics and arranging the Berry modified Bloch equations, Eq. ©, to this form. In section 3 
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we will present the divergence analytically to show that the system is compressible. In section 4 we evaluate 
the Jacobi identity to show that a local canonical form exists and in section 5 we reveal the evolution of 
a probability distribution developed directly from the general form of Darwinian dynamics. Section 6 is a 
summary and discussion. 

2 Evolutionary Decomposition and Conservation of "Energy" 

2.1 Darwinian dynamics 

The Darwinian dynamics arises from a generic nonequilibrium process common in biological, physical, 
and social sciences One of it's main features is to treat all dynamical variables on an equal footing. 
This is typically achieved through expression of the dynamics in a set of first order stochastic differential 
equations. 

Given a generic first order dynamic system of states, x, separated into the deterministic, f(x), and 
stochastic, C(x, t), components: 

x = f(x)+C(x,i) (6) 
The noise is typically approximated as Gaussian and white with zero mean, (C( x > t)) = 0, and variance as 

(C(x,t)C(x,0) = 2o; J D(x) ( 5(t-0. (7) 

Here the nonnegative constant u plays the role of temperature in physics. Further characterization of the 
noise comes from the diffusion matrix D. 

There then exists a unique decomposition as follows [91 I1UI ITT1 112j : 

[S(x)+T(x)]x = -V x $(x)+£(x,i) (8) 

Here S is a symmetric and positive definite friction matrix, related to the zero-mean Gaussian and white 
noise as follows 

(a^t)C(^t')) = 2u;D( X )5(t-t'). (9) 

T(x) is an antisymmetric matrix, ^(x) a scalar potential function and S a symmetric matrix. By definition 
S is positive semi-definite, x r Sx>0, and will have a dissipative effect on the potential function, showing a 
tendency to approach the potential minima. The T matrix describes a non-dissipative part of the dynamics, 
x T Tx = 0, and tends to conserve the potential function. 

2.2 Reformulation of Bloch-Peierls-Berry equations 

Now we reformulate the Bloch-Peierls-Berry equations from the point of view of Darwinian dynamics, 
treating position and momentum coordinates in equal footing. We begin the decomposition of these 
equations (j3J) by arranging them in terms of the complete state vector x in the form of Darwinian dynamics, 
Eq.(jHl). The friction matrix S may not be known. We assume its existence and may take it to be zero at 
the end of calculation whenever needed. The existence of electron-phonon interaction and other dissipative 
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processes and their small possible effect in a solid justifies such a procedure. Beginning from this formation 
of our descriptive matrices a potential function will then become apparent. 
We rewrite the Bloch-Peierls-Berry equations, Eq.©: 



x 



lV k e(k) 



M(x)x 



(10) 



Here the matrix M £ K 6x6 contains the effects of both magnetic field B in position space and Berry 
curvature fi in momentum space. 
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(11) 



With 7 as an identity matrix, Eq.^UJ becomes: 

(/ + M(x))i 




(12) 



This is almost in the form of Eq.(jHJ) of Darwinian dynamics and the right hand side appears similar to 
the gradient of the Bloch Hamiltonian H, Eq. ([T]). Suggesting that the Hamiltonian or energy function 
in the Bloch-Peierls-Berry dynamics may indeed be the original Hamiltonian of Eqs. ©• This makes a 
notable potential function because it is a straightforward representation of the total energy in the system, 
we neglect the constants e and h when speaking of total energy because they can simply be absorbed into 
x. 

We also must mind the order and sign of the state vector x compared to the right side of Eq. (|12|). To 
be consistent an orthogonal transformation matrix R is applied to both sides of the equation (|12|) • 
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The right side of Ea.([12|) becomes 



R 



^V k 6(k) 

"|V r 0(r) / 



Uv k e(k) / 
V x H(x) 



(14) 



4 



and the left side contains an antisymmetric matrix T: 

i?[I + M(x)]x = T(x)x 

where 



(15) 
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Collecting all terms the Bloch-Peierls-Berry equation, Eq.(J3J), is then transformed into the form of Dar- 
winian dynamics, Eq.©: 

T(x) x = -V x W(x) (17) 

The S matrix, as well as the diffusion matrix D, is null in this case, which may be thought as the zero 
electron-phonon interaction limit. The equivalent form to Eq.lJHJ of the Darwinian dynamics is 



x = -Q(x)V x H(x) , 
with Q(x) = T _1 (x). A Poisson bracket can be easily defined as 

df dg 



[/(x), 5 (x)]= QiA< 
»J=i 



dxi dxj 



(18) 



(19) 



Here / and g are two arbitrary functions of the phase space x. The Poisson bracket is anti-symmetric 
because Q is. With this Poisson bracket, the Bloch-Peierls-Berry equations are now 



x = [ft(x),x] , 



(20) 



the familiar form in Hamiltonian dynamics. 



2.3 Conservation of energy 

The potential function of Bloch-Peierls-Berry dynamics in the form of Darwinian dynamics may be clearly 
identified as the Bloch Hamiltonian: 



$(x)=W(x) =-[e(k)+e<Kr)] ■ 



(21) 



It is straightforward to verify that such an "energy" given by the Bloch Hamiltonian is conserved in the 
Bloch-Peierls-Berry dynamics, as expected: 



d_~,{ \ - gx dU{yi) 
dt ~ dt dx 

= x T V x W(x) 

= x T T(x)x 

= 



(22) 
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Thus, we may indeed identify the Bloch Hamiltonian given in Eq.(2) as the candidate for the Hamiltonian 
for the Bloch- Peierls-Berry dynamics of Eq.(|3j). upon the clarification of two questions raised in the liter- 
ature, the validity of the Liouville theorem, and the validity of the Darboux theorem, to be discussed in 
the next two sections. 



3 Compressibility and No Liouville Theorem 

Here we will determine an analytical result for the divergence of the Bloch- Peierls-Berry equations (J3J) from 
the Darwinian dynamics given by Eq. (|17j) . and show that it is indeed non-zero and thus compressible. This 
compressibility feature is pronounced and surprising, as pointed out in Ref.jH]. An explicit demonstration 
is given in light of the present Darwinian dynamics formulation. 
The equation in question is, following Eq.JHJ): 

V x -x = -V x -(T' 1 (x)V x H(x)) (23) 

We first note that because T _1 is antisymmetric only the divergence of the T _1 matrix on the right hand 
side of the equation need be considered: 

V x -x = (V x -T~ 1 )V x W(x) (24) 

Next, we note a useful identity expressing the divergence of T" 1 as a function of the divergence of T, 
because T is easier to operate on. This identity is found as follows, where / is the identity matrix: 

V • / = V x • (T -1 T) 

6 f)T 

= (v x .r- 1 )r + ^r" 1 (^0^- (25) 

= 

The T~ 1 (z, :) denotes the ith row in the T" 1 matrix. 
Solving for V x • T" 1 in Eq. (jlij) we obtain: 



v x .T-^-(Er-(,o£)r-> 

x i=l / 



(26) 



This simplification has also allowed us to express the divergence independent of Ti once combined with 
Eq. (|24|) as follows (this will prove to be useful later.) 



V x -x 



T < 27 > 
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After fully expanding and condensing the operations in parentheses above we obtain a useful and 
simplified expression: 
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(28a) 



(28b) 



(1 + fB 

Eq. (|28b|) is revealed from the fact that the divergence of the Berry curvature, f2, as well as the magnetic 
field, B, are equal to zero. Both are defined as the curl of a vector function and the divergence of a curl 
must be zero, V • V x f = 0: 



V r • B 



V k • SI 



dBi dB 2 dB 3 



dr 2 



<9r 3 

dn 3 







dk± dk 2 dk 3 



(29a) 
(29b) 



In the first row of (|28a|) . for example, — ^c^-fii — Wr-^i is equal to simply jE 1 ^!- Using the solution in 



28b|) our phase divergence, Eq. (|27|l . reduces to a final condensed form as follows. Here J(x) = 1 + IB • f2. 



Vx . i = _VxJ(x) i (3Qa) 

J (X.) 

= ~InJ(x) (30b) 

d ln(l + |B-fi) (30c) 



eft V ft 

We see that the divergence is time- varying given that both B and ft are non-zero and non-orthogonal. 
Thus our system phase space is, in general, compressible, that is the Liouville theorem does not hold. This 
feature was used in Ref.JS] as the indication that the dynamics were non-Hamiltonian. However, different 
opinions have been voiced [U 03 EE EJ. A related point will be explicitly exposed in the next section. 



4 Jacobi Identity 

In this section we follow a standard procedure in mathematical physics |13| I14j to determine when a system 
can be (locally) regarded as a Hamiltonian system. In this regard, there are two important conditions 
the system has to satisfy. The first is the existence of a Poisson bracket which must be explicitly anti- 
symmetric. The second is on the validity of the Darboux theorem. The incompressible condition, or the 
Liouville theorem, does not appear to be essential in this regard. 

The existence of an antisymmetric Poisson bracket has already been defined in the previous section, 
Ea. (|19|) . which endows the system with a specific dynamical structure which may not be a (local) canonical 
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Hamiltonian system. For the Darboux theorem, which guarantees that a non-degenerate system can be 
expressed in local canonical coordinates, the crucial condition to it's validation is the presence of the Jacobi 
identity. 

We express the Jacobi identity as the sum of the cyclic permutations of the double Poisson bracket. 
Here f(i]), g(i]), and ^,(77) are any functions with continuous second derivatives and 77 is the system space 
in consideration. The Jacobi identity is: 

[f,[g,h}] + [h,[f,g}] + [g,[h,f]}=0 (31) 

Direct evaluation to verify the Jacobi identity would involve a great deal of algebra. By computing the 
Poisson bracket as a summation over each vector field we notice most of the terms cancel out. 

i,j=0 3 

Here T~ l {i,j) = Qjj(x) is the ijth element of the T~ l matrix. Also to more easily display derivatives a 
notation will be used in which djjdxi is simply fa. Evaluating the first double Poisson bracket produces, 

[f,[g,h}}= fiT-HiJ)(9kT-\k,l)hi) . 

i,j,k,l=0 1 , , 

(33) 

= Yl fiT-\iJ)[g k jT-\k^)h l + g k T-\kJ)h ld + g k Tf\k,l)h^ 

i,j,k,l=0 

The next two double Poisson brackets give similar results and because T _1 is antisymmetric the first 
2 terms in the parentheses of each will cancel out |15j . In the simple case in which the transformation 
matrix is antisymmetric and independent of x this will alone satisfy the identity. However, difficulty arises 
because dT -1 /dx{ 7^ and we are left with the third term of each permutation containing the derivative 
of the T -1 matrix as follows: 

n 

[f,[g,h}} + [h,[f,g}} + [g,[h,f]} = f l g k hi{T~\iJ)T- 1 (k,l) + T~ 1 (l,j)Tr l (i,k)+T~ 1 (k,j)Tf 1 (l, 

i,j,k,l=0 

(34) 

Eq. (|3*5T) is arranged to collect the arbitrary /, g, h functions. Because each figuhi is unique the ex- 
pression in parentheses must equal zero for every combination of i,k,l in order to satisfy the identity. 

6 

^T-Hi,j)Tr 1 (A:,/)+r- 1 (Z,i)rr 1 (i,fc)+r- 1 (A : ,j)rr 1 (Z,i) =0 for alii , A:, Z (35) 

Because i,k,l = 1,2,..., 6 there are 6 3 equations that we must prove are equal to zero. However, we 
have been able to express all of them in four general forms through index notation. To show how this is 
possible it is helpful to first state the T~ 1 matrix in index notation. Here each quadrant is a 3 x 3 matrix 
with a denoting row and b column 

r -1 fx) = 1 1 £afec ^ c s ab + ^,B a n b \ 

v/dit(r) I -S ab - ^l a B b -e ahc \B c J ' 



s 



det(T(x)) = ( 1 + |B • O ) . (37) 



with the determinant 

rWTfVll = 11 4- 

When solving Eq. (|35j) with the index notation given in (|36|l a distinct difference in solutions is noticed 
when i, k, I have values of 1, 2, or 3 and when they have values of 4, 5, or 6. To take advantage of this we 
will now improve our notation of each subset of i, k, I. A value between 1 and 3 will be denoted by r and 
between 4 and 6 by k. Particular values in r and k can now be described by integers between 1 and 3 for 
the sake of index notation. In the following solutions these integers are given, again, by a, b, and c which 
are unrelated to those in (|36|) . The specific process to arrive at each solution is straightforward but lengthy 
and is excluded here. The final simplified solution for each subset is as follows. 

(Vj.r^j^fVfc.f!) (38a) 

(r ., h , y = ""ffffiV , ' n) + ^aHytV, • B) - (V, ■ B) (38b) 

(t„,t t ,« !o ) = ^i(V r -B) (38c) 

(WJ - ^g^V, ■ b) + s^ )( v k ■ i» - J^jv, • n) (38d) 



As previously mentioned in Eq. (|29j) the divergence of both Berry curvature in momentum space and 
the magnetic field in position are zero and each subset above will also be zero. It should also be noticed that 
the above equations only accommodate one half of the 6 3 total equations. The rest come from permutations 
in (r, k, k) and (k, r, r) and are similar expressions. 

Thus the Jacobi identity is satisfied for the Bloch-Peierls-Berry dynamics. Hence the Darboux theorem 
holds. Therefore, these dynamics can be locally mapped into a canonical form of Hamiltonian dynamics, 
and may be named Hamiltonian dynamics. In fact a global transformation has even been suggested [3]. 
The reformulation of Bloch-Peierls-Berry dynamics into Darwinian dynamics makes this point evident from 
a completely different perspective. 



5 Fokker-Planck Equation and Equilibrium Distribution 



The general form of the present Darwinian dynamics described by Eq. (JSJ) has been shown to correspond 
to a Fokker-Planck equation describing probability evolution in phase space |3 111! I12j : 

dp(x, t) 



Of 



VI 



[£>(x) + Q(x)] [u;V x + V x $(x)] p(x, t) 



(39) 



Here w is a non-negative constant equivalent to temperature. D is a symmetric matrix and Q is an 
antisymmetric matrix so that 

(d + q) = (s + ry 1 . 

Applying Ea. ()39|) to the Bloch-Peierls-Berry equation as given in Eq. (|17|) we find that D = and 
Q = T , thus there is only an anti-symmetric part and the diffusion matrix D is in the zero limit: 

dp(x, t) 



dt 



Q(x)[wV x + V x $(x)] p(x,t) 



(40) 
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The term on the far right side of the equation can be seen to represent the probability flow divergence, 
V x • (T _1 V x <I>p) =>■ V x -(x/j), by substituting Eq. (|17[). Thus in the trivial case where u> = Eq. 1)40(1 
becomes the standard form of the continuity equation. The term V x • (wV x • T~ 1 )p in Eq. ()40|) represents 
the additional effect on distribution at non-zero temperatures. 

We now provide a plausible argument for another choice of equilibrium distribution other than that 
in Ref.[3|. First, we note that in the anti-symmetric matrix T in Eq. 1)17(1 . the geometric phases due to 
the magnetic field B )161 117j and the novel Berry curvature fl are in equal footing. They all contribute 
in the same manner to the geometric phase for a possible close trajectory in the phase space. A unifying 
description of the dynamics is hence achieved in this Darwinian dynamics formulation. 

Second, in the case that dynamics are dissipative in phase space, for example if there are indeed electron- 
phonon and electron-impurity interactions, then the friction or resistance matrix S, and hence the diffusion 
matrix D, is not zero. Such a dissipative dynamics is very likely to be described by Eq.©. Further, such 
a dissipative dynamics in the presence of Berry phase has been well characterized in condensed matter 
physics, such as the dynamics of topological singularities in superconductors and superfluids JH]- Thus, 
extending Bloch-Peierls-Berry dynamics to include dissipation as in the form of Eg. ((17)) or (jHJ) indeed 
appears to be a natural choice. 

Based on the above considerations we may adopt Ea.()39() as the equation for dynamics of non- interacting 
electrons in the presence of dissipation. The equilibrium distribution for such dynamics immediately reads: 



Pe<?(x) = -^exp j- 



The partition function summing up all probability in phase space may be defined as 



(41) 



u 

The time dependent free energy for such an open system may be defined as 



d « x exp <j } . (42) 



(43) 



m = -u J d 6 x p(x, t) In (^y) + ?e q , 



which always decreases towards the equilibrium value T eq = —u\n.Z. 

We note that if treating the equilibrium distribution p eq as a single particle distribution, one immedi- 
ately notices a difference between the present result and the one obtained in Ref.jHj (c.f. their reference 
[21]) in the same limit. 



6 Discussions 



We have investigated the Bloch-Peierls-Berry dynamics from a nonequilibrium dynamical point of view. 
The Bloch-Peierls-Berry equations can be reformulated into a simple and generic form. Using this re- 
formulation we have explicitly demonstrated the compressibility of Bloch-Peierls-Berry dynamics and the 
embedded Jacobi identity, both were pointed out previously in literature. From the point of view of Hamil- 
tonian dynamics, we explicitly showed that the violation of the "naive" Liouville theorem is not essential. 
At the same time we reached a distribution function similar to what is implied in Ref.[3|. 
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There are two new features in our study that we would like to point out. The first one is that our study- 
is in the classical domain, and is for one electron. No many-body effect, as strongly suggested in Ref.jB], 
has been considered here. Hence, from the condensed matter physics point of view, our work may be 
relevant in the dilute and non-interacting limit, with a relatively high temperature. In this sense, it would 
be further interesting theoretically to see how the Boltzmann equation implicitly discussed in Ref.jS] would 
relate to the present Fokker-Planck equation. In particular, how the electron-phonon type interaction with 
explicit energy dissipation would be incorporated into formulation in Ref.[3]. 

The second feature is that due to the open system nature our probability dynamics formulation is 
necessarily in the domain of the canonical ensemble, where the Boltzmann-Gibbs distribution is emphasized. 
This also implies that there is a preferred and natural choice of phase space. Instead, the discussion of 
Hamiltonian flow in Ref's.[3J |5] appears within the micro-canonical ensemble. It is well known that the 
transition from micro-ensemble to canonical ensemble is not unique. This may be the reason that within 
the present nonequilibrium formulation there is no compelling reason to emphasize the Liouville theorem, 
even though the Hamiltonian structure is evident in both approaches. 

In conclusion, the present exploration from the Darwinian dynamics perspective has further revealed 
the richness of Bloch-Peierls-Berry dynamics. 

We thank D. Xiao and Q. Niu for bringing Ref 's. [1J |2j to our attention. This work was partially 
supported by NIH grant #HG002894. 
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